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(57) Abstract: A calibration factor for a sensor includ- 
ing a material-containing conduit and a plurality of mo- 
tion transducers operative to generate motion signals 
representing motion of the conduit is determined from 
a normal modal dynamic characterization of motion of 
the conduit. The calibration factor may be generated by 
generating a solution to a modal differential equation 
of motion in a modal domain including a plurality of 
normal modes of motion of the conduit, and by gener- 
ating the calibration factor from this solution. The nor- 
mal modal dynamic characterization may characterize 
motion of the conduit as a function of a modal Cori- 
olis term that describes coupling among a plurality of 
normal modes of the conduit responsive to mass flow 
in the conduiL The modal Coriolis term may be deter- 
mined from a mode shape function that describes modal 
motion of the conduit as a function of location on the 
conduit. The invention may be embodied as methods, 
apparatus and computer program products. 
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CALIBRATION OF A CORIOLIS MASS FLOW METER USING NORMAL MODAL ANALYSIS 

5 

Field of the Invention 

The present invention relates to sensors and related methods and computer 
program products, and more particularly, to mass flow measurement methods, 
apparatus, computer program products. 

10 Background of the Invention 

Coriolis sensors typically operate by detecting motion of a vibrating conduit 
that contains a flowing material. Properties associated with the material in the 
conduit, such as mass flow, density and the like, may be determined by processing 
signals from motion transducers associated with the conduit, as the vibrational 

15 modes of the vibrating material-filled system generally are affected by the combined 
mass, stiffness and damping characteristics of the containing conduit and the 
material contained therein. 

A typical Coriolis mass flowmeter includes one or more conduits that are 
connected inline in a pipeline or other transport system and convey material, e.g., 

20 fluids, slurries and the like, in the system. Each conduit may be viewed as having a 
set of natural vibrational modes including, for example, simple bending, torsional, 
radial and coupled modes. In a typical Coriolis mass flow measurement 
application, a conduit is excited at resonance in one of its natural vibrational modes 
as a material flows through the conduit, and motion of the conduit is measured at 

25 points spaced along the conduit. Excitation is typically provided by an actuator, 

e.g., an electromechanical device, such as a voice coil-type driver, that perturbs the 
conduit in a periodic fashion. Mass flow rate may be determined by measuring time 
or phase differences between motion at the transducer locations. Exemplary 
Coriolis mass flowmeters are described in U.S. Patent Nos. 4,109,524 to Smith, 

30 4,491 ,025 to Smith et aL, and Re. 31 ,450 to Smith. 

The accuracy of Coriolis mass flowmeters may be compromised by the 
mountings which constrain vibration of the conduit. The affects of these constraints 
may be reduced by using flowmeter designs that are balanced to reduce effects 
attributable to external vibration, and by using frequency domain filters, e.g., 
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bandpass filters designed to filter out components of the motion signals away from 
the excitation frequency. However, mechanical filtering approaches are often 
limited by mechanical considerations, e.g., material limitations, mounting 
constraints, weight limitations, size limitations and the like, and frequency domain 
5 filtering may be ineffective at removing unwanted vibrational contributions near the 
excitation frequency. 

Conventional Coriolis mass flowmeters typically use a "calibration factor" to 
scale mass flow rate to time difference or phase measurements to generate a mass 
flow estimate. Typically, in order to generate a calibration factor, a calibration or 
10 "proving" procedure is performed in which a material (e.g., water) is passed through 
the vibrating conduit of the mass flowmeter in a test fixture while signals produced 
by the motion transducers of the mass flowmeter are processed to determine time 
or phase differences . This data is then processed to generate an estimated, 
calibration factor for the meter. 

15 This approach can have several disadvantages. The calibration procedure 

can be time-consuming, labor-intensive and costly. In addition, because the: 
calibration is typically not performed in the field, it may be subject to inaccuracies 
due to differences between the test conditions and the conditions of the field, 
including differences in mounting conditions. 

20 Summary of the Invention 

According to embodiments of the present invention, a calibration factor for a 
parameter sensor including a conduit configured to contain a material and a 
plurality of motion transducers operative to generate motion signals representing 
motion of the conduit is determined from a normal modal dynamic characterization 

25 of motion of the conduit. The calibration factor may, for example, relate mass flow 
rate to a spatio-temporal relationship among motion signals produced by the motion 
transducers, such as a time difference or phase relationship. The calibration factor 
may be generated, for example, by generating a solution to a modal domain 
differential equation of motion, and by generating the calibration factor from this 

30 solution. 

According to some embodiments of the present invention, the normal modal 

dynamic characterization characterizes motion of the conduit as a function of a 

modal Coriolis term that describes coupling among a plurality of normal modes of 

the conduit responsive to mass flow in the conduit. Determination of the calibration 

2 
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factor may be preceded by determining the modal Coriolis term from a mode shape 
function that describes motion of the conduit in a normal mode as a function of 
location on the conduit. Modal mass and stiffness terms of the normal modal 
characterization may be generated using, for example, conventional modal analysis 
5 techniques. . 

For example, in some embodiments, the modal Coriolis term may be 
determined from a spatial Coriolis characterization that describes motion of a 
plurality of discrete locations of the conduit in a spatial domain in response to a 
predetermined mass flow. The modal Coriolis term may be determined from the 

10 spatial Coriolis characterization, (e.g., a spatial Coriolis matrix), using a 

transformation that relates the spatial domain to the plurality of normal modes. The 
plurality of discrete locations associated with the spatial Coriolis characterization 
may include locations other than the transducer locations. 

The transformation that relates the spatial domain to the plurality of normal 

15 modes may be generated from a mode shape function that describes motion of the 
conduit in a normal mode as a function of location on the conduit. The mode shape 
function may be determined from a predetermined eigenvalue and a predetermined 
boundary condition, for example, from an assumption as to a boundary condition to 
which the conduit will be constrained. 

20 According to other embodiments of the present invention, an orthogonality of 

a mode shape function for a predetermined mass flow in the conduit is determined 
in order to determine a modal Coriolis term of a normal modal dynamic 
characterization of conduit motion. Similar orthogonality determinations can be 
used to generate modal mass and stiffness terms of the normal modal 

25 characterization. Such terms may also be generated using conventional modal 
analysis techniques. 

In other embodiments of the present invention, an estimated spatial 
response for a plurality of locations of the conduit is generated from a normal modal 
dynamic characterization, and a calibration factor is generated from the estimated 

30 spatial response. For example, a spatio-temporal relationship among movements 
at a plurality of locations of the conduit, such as a time difference or a phase 
relationship, may be determined, and the calibration factor may be determined from 
the spatio-temporal relationship. 
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In still other embodiments of the present invention, mass flow of a material in 
a conduit may be determined by determining a calibration factor from a normal 
modal dynamic characterization of motion of the conduit, generating a plurality of 
motion signals representing motion of the conduit at a plurality of locations on the 
5 conduit, and processing the motion signals according to the determined calibration 
factor to produce a mass flow estimate. For example, the calibration factor may be 
determined from a representation of a modal differential equation of motion, as 
described above. The mass flow estimate may be generated, for example, by 
determining a spatio-temporal relationship among the plurality of motion signals 

10 and applying the calibration factor to the determined spatio-temporal relationship to 
generate a mass flow estimate. 

According to aspects of the present invention, the calibration factor may be 
determined from a mode shape function that describes motion of the conduit in a 
normal mode as a function of location on the conduit. The mode shape function 

15 may be based on an assumed boundary condition for the conduit. Motion of the 
conduit may be constrained to induce a boundary condition that approximates the 
assumed boundary condition, for example, by attaching the conduit to structure or 
by applying a force to the conduit responsive to motion of the conduit to induce a 
boundary condition that approximates the assumed boundary condition. 

20 In still other embodiments of the present invention, an apparatus for 

calibrating a parameter sensor including a conduit configured to contain a material, 
a plurality of motion transducers operative to generate motion signals representing 
motion of the conduit and a mass flow estimator circuit operative to estimate mass 
flow from the motion signals according to a calibration factor is provided. The 

25 apparatus includes a calibration factor generator circuit operative to generate a 
calibration factor from a normal modal dynamic characterization of motion of the 
conduit, the calibration factor generator circuit including an interface circuit 
operative to convey the generated calibration factor to the mass flow estimator 
circuit. 

30 In some embodiments, the calibration factor generator circuit includes a 

normal modal dynamic characterizer circuit operative to generate a normal modal 

dynamic characterization of motion of the conduit and a calibration factor 

determiner circuit operative to determine the calibration factor from the normal 

modal characterization. The normal modal characterization may be, for example, a 

4 
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representation of a modal differential equation of motion in a modal domain 
including a plurality of normal modes of motion of the conduit, and the calibration 
factor determiner circuit may be operative to generate a solution of the modal 
differential equation of motion from the representation of the modal differential 
5 equation of motion and to generate the calibration factor from the solution of the 
modal differential equation of motion using, for example, the techniques described 
above. The calibration factor generator circuit may be integrated with the mass 
flow estimator circuit in, for example, a mass flow sensor apparatus, or may be 
separately implemented in, for example, test or process control equipment. 

10 In yet other embodiments of the present invention, a mass flow sensor 

includes a conduit, such as a substantially straight tube, that is configured to 
contain a material. The mass flow sensor further includes an actuator operative to 
excite the conduit, a plurality of motion transducers operative to generate motion 
signals representing motion of the conduit and a mass flow estimator circuit 

1 5 configured to receive the motion signals and operative to generate a mass flow 
estimate therefrom using a calibration factor derived from a normal modal 
characterization of motion of the conduit that assumes a predetermined boundary 
condition for the conduit. Means are provided for constraining motion of the conduit 
to approximate the predetermined boundary condition. 

20 The means for constraining conduit motion may include, for example, means 

for attaching the conduit to a structure, such as clamps, welds or other fastening 
means that attach spaced apart locations of the conduit to a rigid structure. The 
means for constraining motion of the conduit may also include means for applying a 
force to the conduit responsive to motion of the conduit to thereby induce a 

25 boundary condition that approximates the assumed boundary condition. For 
example, the means for applying a force to the conduit may include at least two 
motion transducers operative to generate motion signals representing motion of the 
conduit, a plurality of actuators operatively associated with the conduit, and a 
shape control circuit configured to receive motion signals from at least two motion 

30 transducers and operative to drive the plurality of actuators responsively thereto. 
According to other embodiments of the present invention, a computer 
program product for characterizing a parameter sensor including a conduit 
configured to contain a material and a plurality of motion transducers operative to 
generate motion signals representing motion of the conduit is provided. The 
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computer program product includes computer-readable program code embodied in 
a computer-readable storage medium, the computer-readable program code 
including program code for determining a calibration factor for processing motion 
signals produced by the plurality of motion transducers from a normal modal 
5 dynamic characterization of motion of the conduit. The program code for 

determining a calibration factor may include program code for generating a solution 
to a modal differential equation of motion in a modal domain including a plurality of 
normal modes of motion of the conduit, and program code for generating the 
calibration factor from the generated solution of the modal differential equation of 
10 motion. The normal modal dynamic characterization may characterize motion of 
the conduit as a function of a modal Coriolis term that describes coupling among a 
plurality of normal modes of the conduit responsive to mass flow in the conduit. 
The computer readable program code may include program code for determining 
the modal Coriolis term from a mode shape function that describes motion of the 
15 conduit in a normal mode as a function of location on the conduit. _ - r. 

Brief Description of the Drawings 
FIG. 1 is a schematic diagram of a mass flow sensor apparatus according to 
embodiments of the present invention. 

FIG. 2 is a schematic diagram of a mass flow sensor apparatus according to 
20 other embodiments of the present invention. 

FIG. 3 is a schematic diagram of a calibration factor generator circuit 
according to embodiments of the present invention. 

FIGs. 4-6 are flowcharts illustrating exemplary operations for determining a 
calibration factor for a mass flow sensor according to embodiments of the present 
25 invention. 

FIGs. 7-11 are schematic diagrams mass flow sensor apparatus according 
to various embodiments of the present invention. 

FIG. 12 is a flowchart illustrating exemplary operations for estimating mass 
flow according to embodiments of the present invention. 
30 FIG. 13 is a free-body diagram illustrating behavior of a material-containing 

conduit. 

Detailed Description of Embodiments 

The present invention now will be described more fully hereinafter with 

reference to the accompanying drawings, in which embodiments of the invention 

6 
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are shown. This invention may, however, be embodied in many different forms and 
should not be construed as limited to the embodiments set forth herein; rather, 
these embodiments are provided so that this disclosure will be thorough and 
complete, and will fully convey the scope of the invention to those skilled in the art. 
5 Like numbers refer to like elements throughout. As will be appreciated by one of 
skill in the art, the present invention may be embodied as systems (apparatus), 
methods, or computer program products. 

Embodiments of the present invention described herein relate to Coriolis 
mass flowmeters that utilize a single, substantially straight conduit that does not 
10 require, for example, a mechanically-tuned balance beam. The present invention 
may be particularly advantageous in such applications. However, those skilled in 
the art will appreciate that the present invention can be used with other types of 
flowmeter structures, including conventional curved tube structures and straight 
tube structures including mechanical balance beams. 
1 5 As will be appreciated by one of skill in the art, the present invention may be 

embodied as an apparatus and/or method and/or computer program product. The 
present invention may be implemented in hardware or in a combination of hardware 
and software aspects. Furthermore, the present invention may also take the form of 
a computer program product including a computer-usable storage medium having 
20 computer-usable program code embodied in the medium. Any suitable computer 
readable medium may be utilized, including semiconductor memory devices (e.g., 
RAMs, ROMs, EEPROMs, and the like), hard disks, CD-ROMs, optical storage 
devices, and magnetic storage devices. 

Computer program code for carrying out operations of the present invention 
25 may be written in an object oriented programming language, such as Java® or 
C++, and/or in a procedural programming languages, such as "C." The program 
code may execute on a single computer or data processing device, such as a 
microcontroller, microprocessor, or digital signal processor (DSP), or may be 
executed on multiple devices, for example, on multiple data processing devices that 
30 communicate via serial or parallel data busses within an electronic circuit board, 
chassis or assembly, or which form part of a data communications network such as 
a local area network (LAN), wide area network (WAN), or internet. 

The present invention is described below with reference to flowchart 

illustrations and/or block diagrams of methods, apparatus (systems) and computer 

7 
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program products according to embodiments of the invention. It will be understood 
that blocks of the flowchart illustrations and/or block diagrams, and combinations of 
blocks in the flowchart illustrations and/or block diagrams, can be implemented by 
computer program code (instructions). These computer program code may be 
5 provided to a processor of a general purpose computer, special purpose computer, 
or other programmable data processing apparatus to produce a machine, such that 
the instructions, which execute via the processor of the computer or other 
programmable data processing apparatus, create means for implementing the 
functions specified in the flowchart and/or block diagram block or blocks. These 

10 computer program products also may be embodied in a computer-readable storage 
medium (e.g., magnetic disk or semiconductor memory, code magnetic memory or 
the like) that can direct a computer or other programmable data processing 
apparatus to function in a particular manner, such that the computer program 
stored in the computer-readable memory produce an article of manufacture 

1 5 including instruction means which implement the function specified in the flowchart 

and/or block diagram block or blocks. The computer program code may also be ^ 
loaded onto a computer or other programmable data processing apparatus to ^ 
cause a series of operational steps to be performed on the computer or other 
programmable apparatus to produce a computer implemented process such that 

20 the code that executes on the computer or other programmable apparatus provides 
steps for implementing the functions specified in the flowchart and/or block diagram 
block or blocks. 
Overview 

According to embodiments of the present invention, a Coriolis mass flow 

25 sensor may be calibrated using a normal modal dynamic characterization of the 

motion of the sensor's material-containing conduit. In some embodiments of the 

present invention, this normal modal dynamic characterization comprises a modal 

differential equation of motion that includes a Coriolis term that represents coupling 

among normal modes of the flow sensor conduit in response to mass flow. 

30 According to some embodiments of the present invention, the Coriolis term 

comprises a modal Coriolis matrix that is generated from a spatial (or "physical") 

Coriolis matrix that describes mass flow dependent motion of the conduit in a 

spatial domain (e.g., Cartesion). The modal Coriolis .matrix may be derived from 

the spatial Coriolis matrix using a transformation that is derived from a mode shape 

8 
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function that describes modal behavior of the conduit as a function of location on 
the conduit. The mode shape function may be generated from motion signals 
generated by motion transducers operatively associated with the sensor conduit In 
particular, discrete modal response data from the transducer locations may be 
used to generate a continuous mode shape function, which may then be used to 
calculate a modal-to-spatial transformation ® virttial that includes "virtual" locations of 

the conduit other than the actual transducer locations, such as locations between 
the actual transducer locations. 

The spatial Coriolis matrix may represent a combination of local Coriolis 
matrices that correspond to beam like segments defined by nodes corresponding to 
the virtual locations. The local Coriolis matrices represent motion of these 
segments under a predetermined mass flow condition, (e.g., a predetermined mass 
flow rate). The transformation <X> Wr/wa/ may be applied to the spatial Coriolis matrix 
to determine a modal Coriolis term (matrix) of a (vector) modal differential equation 
of motion that also includes modal mass, modal damping and modal stiffness 
terms. 

The modal differential equation of motion may be solved io determine a 
modal response at the predetermined mass flow rate from a known input/excitation. 
The modal response may be converted to a spatial domain response at the actual 
transducer locations, and a calibration factor may then be determined based on the 
estimated spatial domain response and the predetermined mass flow rate. This 
calibration factor may then be applied to motion signals generated in response to 
an unknown mass flow to generate an estimate of the unknown mass flow. 

According to other embodiments of the present invention, a more 
generalized approach to determining a modal Coriolis matrix may be taken. In 
particular, instead of determining a modal Coriolis matrix from a discrete 
transformation of a spatial Coriolis matrix, a modal Coriolis matrix may be 
generated by determining an orthogonality of a mode shape function at a 
predetermined mass flow rate. A modal Coriolis matrix so determined may be then 
used to generate a solution to a modal differential equation of motion for the 
conduit and, from this, a calibration factor may be determined along the lines 
described above. 
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Embodiments of the present invention offer several potential advantages 
over conventional calibration techniques. In particular, use of calibration 
techniques according to embodiments of the present invention may eliminate a 
need to pass material through the mass flow sensor, as is commonly required in 
conventional proving processes. This can provide, for example, an in situ 
calibration capability. 

The present invention can also facilitate the practical implementation of 

straight tube mass flow sensors such as "clamp-on" mass flow sensors, e.g., 

sensors that include motion transducers, actuators and related components that 

are configured to operatively engage to a pipeline such that the pipeline itself 

serves as part of the sensor. A potential difficulty in calibrating such a clamp on 

sensor is that, unlike a conventional Coriolis mass flow meter, the flow conduit may 

not be available to the sensor's manufacturer. Although one theoretically could 

connect a prover to the pipeline and calibrate a clamp-on flow sensor against a 

master meter, such a procedure may be prohibitively time consuming, costly, or 

impractical due to difficulty in accessing the sensor, harshness of the pipeline 

environment, and the like. According to embodiments of the present invention, in 

situ calibration of such a sensor can be performed by measuring the modal 

parameters of the installation after the sensor components are installed. This can 

eliminate the need for a calibration facility or master meter. 

A potential advantage of such a clamp-on flow sensor is that it may be 

moved to different locations and used, for example, as a temporary sensor for 

trouble shooting purposes. As each pipeline the sensor components are attached 

to may differ, a calibration may be needed each time the sensor is relocated. In 

situ calibration according to embodiments of the invention can provide a practical 

and cost-effective way to perform such a calibration. 

FIG. 1 illustrates a mass flow sensor apparatus 100 according to 

embodiments of the present invention. The sensor 100 includes a conduit 103 

configured to contain a material 108. An actuator 106 is operative to excite the 

conduit 103. For example, the actuator 106 may comprise an inertial or relative 

actuator driven by a drive circuit (not shown) that causes the actuator 106 to vibrate 

the conduit 103 at a predetermined frequency; however, it will be appreciated that 

other mechanisms for imparting vibratory motion to the conduit 103 may be used. 

Motion transducers 105 (e.g., inertial velocity transducers, accelerometers or other 

10 
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motion-sensing devices) are positioned along the conduit 103 and produce motion 
signals 107 representing motion of the conduit 103 in response to a plurality of 
forces that may include, for example, a drive force F D imparted by the actuator 106, 
Coriolis forces F c arising from the flowing material 108, and extraneous forces Fe. 
5 The sensor apparatus 100 further includes a mass flow determiner circuit 

110 that receives the motion signals 107 generated by the motion transducers 105. 
The mass flow determiner circuit 1 10 is operative to process the motion signals 107 
according to a calibration factor 155 generated by a calibration factor generator 
circuit 150. The calibration factor generator circuit 150 determines the calibration 

10 factor 155 from a normal modal dynamic characterization of motion of the conduit 
103, for example, a modal differential equation of motion as described in detail 
below. As also described in detail below, the normal modal dynamic 
characterization may be derived from a mode shape function generated from 
discrete modal response data, e.g., response data generated by motion signals 

15 representing motion of the conduit, such as motion signals generated by the motion 
transducers 105. 

FIG. 2 illustrates a mass flow sensor apparatus 200 according to other 
embodiments of the present invention. The sensor 200 includes a conduit 203 
configured to contain a material 208 A relative actuator 206 is attached to a 

20 structure 204 (which may comprise, for example, a housing that encloses the 

conduit 203) such that it induces motion of the conduit 203 relative to the structure 
204. As shown, relative motion transducers 205 (e.g., velocity transducers, 
accelerometers or other motion-sensing devices) are positioned along the conduit 
203 and measure motion of the conduit 203 relative to the structure 204. The 

25 motion transducers 205 produce motion signals 207 representing motion of the 

conduit 203 in response to forces exerted on the conduit 203. !t will be understood 
that, although the motion transducers 205 are shown as relative motion 
transducers, absolute or "inertial" motion sensors (e.g., accelerometers), or a 
combination of relative and absolute motion transducers (e.g., a combination of 

30 relative velocity transducers and accelerometers), may be used instead. 

The sensor apparatus 200 further includes a mass flow determiner circuit 
210 that receives motion signals generated by the motion transducers 205. The 
mass flow determiner circuit 210 is operative to process the motion signals 
according to a calibration factor 255 generated by a calibration factor generator 
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circuit 250. The calibration factor generator circuit 250 determines the calibration 
factor from a normal modal dynamic characterization of motion of the conduit 203. 

It will be understood that, in general, functions of a calibration factor 
generator circuit according to embodiments of the present invention, such as the 
calibration factor generator circuits 150, 250 of FIGs. 1 and 2, may be implemented 
in analog circuitry, digital circuitry and combinations thereof. This circuitry may be 
integrated in a single device, such as a microprocessor, microcontroller, digital 
signal processor (DSP), application-specific integrated circuit or other monolithic 
device, or may be distributed across multiple devices. It will be further appreciated 
that a calibration factor generator circuit according to embodiments of the present 
invention may be integrated within a flow sensor apparatus, e.g., using processing 
circuitry also used to implement a mass flow determiner circuit such as the mass 
flow determiner circuits 1 10, 210 of FIG. 1, or may implemented in separate 
hardware, such as a piece of test or maintenance equipment or process control 
computing apparatus. For example, a calibration factor generator circuit could.be 
implemented in test equipment designed to interface with motion transducers of a 
flow sensor, to determine a calibration factor from motion signals generated by the 
motion transducers, and to download the determined calibration factor to a mass 
flow determiner circuit of the flow sensor. 

The following discussion provides an underlying theoretical description that 
supports apparatus, methods and computer program products according to 
embodiments of the present invention. It will be appreciated that the scope of the 
present invention is not limited to this theoretical explanation. For example, the 
discussion that follows describes computational operations. Although these 
operations are presented with a particular sequence and organization to facilitate 
understanding of the present invention, it will be appreciated that the computational 
operations may be reordered, combined in a different fashion or otherwise modified 
within the scope of the present invention. Consequently, it will be further 
understood that, in general, the present invention includes not only the specific 
computational operations described herein, but also equivalent operations. 

A modal differential equation of motion for a single straight flowmeter 
conduit, such as the conduit 103 of FIG. 1, is: 

12 
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[mr] {//} + [cr] {77} + [dr] {77} + [kr] {/?} = {N} , 



(D 



where {7}, {77}, {77} represent displacement, velocity and acceleration vectors, 
respectively, in a normal modal domain comprising a plurality of normal (i.e., 
5 orthogonal) modes of vibration of the conduit, [mr] is a modal mass matrix, [cr] is a 
modal Coriolis matrix that describes coupling between the normal modes of the 
conduit responsive to mass flow in the conduit, [dr] is a modal damping matrix, [kr] 
is a modal stiffness matrix, and {N} is a modal excitation applied to the conduit. 
The modal mass matrix [mr] and the modal stiffness matrix [kr] may be determined 
10 using conventional modal analysis techniques. However, determining the modal 
Coriolis matrix [cr] can be problematic. As described in greater detail below, a flow 
calibration factor (FCF) may be assumed to be relatively insensitive to damping and 
strongly correlated to the modal Coriolis matrix [cr]. 



15 "DETERMINING PROPERTIES OF A FLOW TUBE AND OF A FLUID FLOWING 
THROUGH A FLOW TUBE OF A CORIOLIS FLOWMETER," filed concurrently 
herewith and incorporated herein by reference In its entirety, describes how 
discrete mode shape information for a flowmeter conduit gained from motion 
transducers that generate motion signals representative of motion of the flowmeter 

20 conduit can be used to determine a continuous mode shape function, i.e., a 

function that describes motion of the conduit in a plurality of normal modes as a 
function of location on the conduit. From such a continuous mode shape function, 
"virtual" pickoff responses, i.e., responses for locations other than the transducer 
locations, may be synthesized. 

25 The aforementioned Patent Application 09/941 ,332 describes continuous 

mode shape functions of the form: 



A United States Patent Application Serial Number 09/941,332, entitled 




bc x 
bc 2 
bc 3 



r 



(2) 
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wherein <D r is a mode shape function for an r-th normal mode, and wherein 
eigenvalues /Ll,.,A2 r and boundary conditions {bci, bc 2) bc 3 , bc 4 } r of the mode 
shape function 4> r may be calculated from modal measurements. 

Such a mode shape function can be used to transform a spatial or "physical" 
Coriolis matrix to a modal Coriolis matrix. In particular, a modal transformation 
[^wrw] for a plurality of selected locations of a flowmeter conduit may be 
synthesized using mode shape functions as described above. The number of such 
"virtual" locations preferably is greater than the number of locations at which motion 
transducers are present. A modal Coriolis matrix [cr] may then be calculated from a 
spatial Coriolis matrix [C] for the plurality of virtual locations using [O^,] 
according to the following: 



[crM® nmol ] T [C)[<S> vimal ] 



(3) 



The spatial Coriolis matrix [C] may be generated using conventional 
techniques. The plurality of virtual locations of the conduit may be viewed as 
defining a plurality of beamlike segments having first and second nodes that 
correspond to virtual locations, with a spatial Coriolis matrix for each segment at a 
predetermined mass flow rate m being given by: 



m 



I0gzl + l 2 
60gz + 5l 2 

-I 

iogzi+r- 



I0gzl + l 2 
60gz + 5l 7 

0 

10gzl + l 2 
~60gz + 5L 2 
l 4 



1 

I0gzl + l 2 
60gz + 5l 2 

0 

I0gzl + l 2 



60gz + 5l 2 360gz + 30/ 2 60gz + 5l 2 



I0gzl + l 2 
60gz+5! 2 
-/ 4 

360gz+30/ 2 
10gzl + l 2 
60gz + 5! 2 

0 



(4) 



wherein the first row/column of the "local" Coriolis matrix [C| 0ca i]corresponds to the 
physical displacement of a first node of a segment, the second row and column 
corresponds to the angular displacement of the first node of the segment, the third 
row/column corresponds to the physical displacement of a second node of the 
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segment, and the forth row and column corresponds to the angular displacement of 
the second node of the segment. In equation (4), / is the length of the segment and 
gz represents a combination material and geometric properties of the conduit, 
which may be given by: 



wherein E is the modulus of elasticity for the conduit material, l t is the area moment 
of inertia for the conduit, G is the shear modulus of elasticity for the conduit 

material, A t is the cross-sectional area of the conduit, k is a numerical factor 
dependent on the cross-sectional shape of the conduit. The spatial Coriolis matrix 
[C] for the conduit may be constructed from the local Coriolis matrices [C| OC ai], using 
techniques known to those skilled in the art. 

A modal dynamic stiffness matrix [MDS] for a drive frequency excitation, /.e., 
a frequency at which the conduit is vibrated, can then be created and inverted to 
determine a moda! frequency response function [MDS]" 1 that can be used to 
determine a modal response {77} for the predetermined mass flow rate m : 

{v} = W + jm^ [cr] + ja> arv [dr] + [Ar])"' W = [MDS]' 1 {N} . (6) 

The modal response {7} can be transformed into a spatial domain response vector 
{y} as follows: 

{y} = [V^IMDS]- 1 [0„,f {F} . (7) 

A spatio-temporal relationship between transducer locations, such as a time 
difference ("delta-t") or a phase relationship between motion signals corresponding 
to the transducer locations, for the predetermined mass flow rate m can be 
determined from the estimated spatial response {y}. For example, a spatial 
integration technique for generating a time difference or phase estimate from 
spatial motion information provided by a plurality of spatially distributed motion 
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transducers may be used. An example of a spatial integration technique is 
described in United States Patent No. 6,233,526, incorporated herein by reference 
in its entirety. 

The time difference or phase relationship, which may be referred to generally 
5 as A, and the predetermined mass flow rate m can be used to determine a flow 
calibration factor FCF in the appropriate units (e.g., mass/sec/sec or 
mass/sec/degrees): 



10 



15 



FCF=—. (8) 
A 



A time difference or phase relationship measurement (generally, A measured ) 
generated responsive to an unknown mass flow rate m unhlown may be multiplied by 
the calibration factor FCF to generate an estimate of the unknown mass flow rate 

^unknown * 



= FCF X ^measured " A o) > (») 



where A 0 is an offset (e.g., time or phase) between transducer signals at 

substantially zero flow ("zero offset"). 

20 Experimentation indicates that, in many practical applications, mass flow rate 

is nearly linearly proportional to time-difference over a range of mass flow rates, 
such that FCF may be assumed to be constant. However, if FCF is not sufficiently 
constant with flow, it could be estimated over a range of assumed flow rates to 
develop, for example, a lookup table or functional representation thereof. FCF 

25 could also be modified based on other parameters, such as temperature, e.g., one 
could compensate for a change in FCF due to changes in modulus (E) with 
temperature. 

If it is assumed that the ends of a flow conduit of a mass flow sensor are 
constrained such that forced response of the motion of the conduit ends is about 
30 zero, and if it is further assumed that the case or other structure supporting the 
motion transducers is stiff and massive enough such that it acts like a rigid body, 
the motion signals generated by the sensor's motion transducers can be treated as 
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approximating absolute velocity signals. Under such conditions, the flow conduit 
can be treated as having approximately fixed boundary conditions at the frequency 
of excitation of the conduit The boundary conditions and eigenvalues may, 
therefore, be determined, at least for Euler beam assumptions. 

For a flow sensor conduit 1310 as illustrated in the free-body diagram of FIG. 

13: 



V 9 - p t A t - dy + Fc + rc(0 + a) +F(xJ) = 0 
M' — Sy' + V- p,I t — Mc = 0 



10 where V is the shear force exerted on the conduit 1310, p t \s the density of the 

material of the conduit 1310, >4/is the cross-sectional area of the conduit 1310, d is 
a damping coefficient, y is the acceleration of the conduit 1310, Fc is force coupling 
between a material 1320 contained in the conduit 1310, ris viscose shear force 
between the material 1320 in the conduit and the inner wall of the conduit 1310, 6 is 

1 5 the angular displacement of the conduit 1310, a is the shear distortion of the 

conduit 1310, F(x, i) is the drive force applied to the conduit 1310, M is the bending 
moment of the conduit 1310, S is the tension acting on the conduit 1310, is the 
inertia per unit length of the conduit 1310, and c is circumference of the inner 
surface of the conduit 1310, and wherein dots indicate time derivatives and primes 

20 indicate spatial derivatives with respect to length. For the material 1320 contained 
within the conduit 1310: 



( d d V 
-p f I f Q +Mc + PA f y' = 0 



25 where pAf is the mass per unit length of the material 1320 in the conduit 1310, y is 
the displacement of the conduit 1310, vy is the velocity of the material 1320 in the 
conduit 1310, pff f \s the inertia per unit length of the material 1320 in the conduit 
1310, 0 is the angular acceleration of the conduit 1310, P is the pressure of the 
material 1320 in the conduit 1310, and >4fis the cross-sectional area of the material 

30 1 320 in the conduit 1310. 
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Eliminating coupling forces (Fc), moments (Mc) and fluid shear (□) on the 
wall of the conduit "beam": 



V -p,A- + P/Af [j f + v, J^) y + F{x, 0 = 0 
M' — Sy' + V- (p,I, + p f I f ) 9 + PA f y' = 0 



(12) 



From basic elasticity principles: 

y'(x,t) = 0(x,t) + a(x,t) V(xj) = kGA,a{x,t) M (x, t) = EI ,9' (x, i) (13) 
1 0 This can be reduced to two equations in terms of y(x,t) and Dx.t. 



kGA t (y - 0') - Pl A,'y ~ dy - Pf A, ^- + v f —J y+ f(x, t) = 0 - -~ 
EI ,9" - Sy' + JcGA, (/ -0)-(p,I l + p f I f ) 9 + PA f y' = 0 



Non-dimensionalizing the basic independent variables in equation (14), 
1 5 where the under-bar indicates a dimensionless variable of space or time and the 
pointed brackets indicate the U.S. Customary units for the variables: 
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x = Lx dx = Ldx 



/ - — I dt = — d± co r =— ^= — 
co x co x — co x co { 



4» r (x) = O r (x) = O £ (x)<0) 
d >r (x) = [ e A1 ' j: e"*''* e* 2 '* e"^ v ]{£C r }(0> 
O r (x) = [e A, ' £j£ e y/t2 ' ias e- JX2 ' Lx ]{BC r }(0) 

<&,(*) =[V' i2 e M ' Lx - e m ' Lx - e- JX2 ' L£ ]{BC r }{0) 

* r (*) = ¥ r (*) = l^(*)(l/m) 

^ r (x)=[e XKx e -u,* e J*2 r * e -^^][ Mr ]{ jB C r }{l//n> 
¥ r (x)=[e M ' L£ e"^^ e y ^ ££ e -^' i£ ][« r ]{5C r }(l/m> 
T,(x) = [e' lI " tJE e^ 1 * e JX2 ' Lx e jX2 ' Lz -\i±\{BC r }(0) 
[«J = l[« r j 

y{x,t) = £<D r (x>/ r (t) ^ 2> r {x)rj r (t) = y(x,t) = (in) 

r r 

y&t) - t2>. UK 00 = (*k U)(o) 

^ r r 

*<x,f) = 2> r (*) 7r (0 = 007r (0 = = (°) 

r r 



Differentiating and integrating for the non-dimensionalized mode shapes ^r,lFJ-: 



(15) 
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9/ 1 rV_/ dt 1 rW dt \/sec/ 

a/ 2 w e* a? \/sec7 



gVM _ g, eg, (*) a?, (|) = ag, (x) e%. (/) , 
a*a* l dx dt 1 a* a? \/ sec / 



££M = l£?VM, rW= ±£2lM I7 ,( ( )(y) 



^.^^.^wAO^ (16) 



_ (x) £VM . ^.» r (x) *2k« (>M 

^2 i r v_y ^ 2 i ,._rv-y fl/ \/sec/ 



3x L dx dx L dx x ' 

£ <D r (x) & = Z0O r (x) rfx = (/») (17> 
I>. (*)<■>. {x)dx = L0 L <S> m (x)®„ (x)dx = (in) 



5 In equation (17), the limits of integration x y L and x y L are typically 0 to 1 or -Va to 

+>£ depending on where the coordinate system has been placed. Equation (18) is 
the same as equation (14), except all the non-dimensional relationships have been 
included: 
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kGA, 



EI, 
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EI, ^ W dt 2 EI, dt_ 



2 Pf A f v f <D l d^jx) d^jt) p,A,v/ 1 ?Q,(x) 1 , x_ 

_ H. Fx dt E%~l~~dx r ~ Vr U) EI,V*- d '- 0 



dx dt 



l^, (l)+ (^£2^, rW+ Mr^M_^ w l, rW -... 



a*' 



9x 



£7, — w a* 2 



(18) 



Multiplying through equation (18) by L 2 and L, respectively: 



kGAL 2 



(\d 2 O r (x) ldV r (x) 



EI 



L dx 2 



L dx } 



dt 



2 Pr A f v f ^ i aa>,(x)a^(0 P/ a /V /l 2 i d 2 ^(x) (t)+ e_ F(x t)=0 



Z, dx dt 



EI, L dx' 



EI. 



.^ w+ (^lM i(i)+ ^(lM_^ (£ ^ (!) -. 



z a* 

/?/<y, 2 Z 2 1 



zz, z— w a/ 2 



(19) 



In equation (19), S is positive for tension, therefore S is negative when the tension 
10 is positive: 



— EI, EI, — EI, 

_ p f A f v/L 2 D _kGAL 2 . PfA f v/ 



EI, 



m = ■ 



EI, " EI, 
EI. 



(20) 



- EI, — EI. V ' — Z:7, 
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Putting equation (19) in compact matrix form: 



' pA 0 " 
0 pi 



U-B) 



s 2 D 5 iD 



d + 2mH 0 



ax 

0 0 



Ml) 



a* 2 

d 



ax 
a 2 D 



v - -'ax ax 2 



(21) 



5 It may be noted that the mode shape functions O r (x) and (x) are non- 
dimensional. The 1/L term is carried along to non-dimensionalize the integration 
when calculating orthogonality. The "box" symbol is a placeholder for the 
differential operator. The first term of equation (21) accounts for all the mass and 
inertia of the beam. The second term accounts for proportional damping and the 

10 Coriolis effects. The third term accounts for the stiffness, centrifugal forces, 
pressure, and tension in the beam. The stiffness matrix is mainly diagonal with 
some slightly skew symmetric terms. 

To calculate the orthogonality, or "self-adjointedness" of any two modes m 
and r , each term of the differential equation may be integrated over the entire 

15 length of the beam: 



r*f®.( x )T 
J *K(*)J 



f*f<">.(*)l 



Tr 



pA 0' 
0 pi 



d+ 2mS. o 
_ — dx 

0 0 



life) 



a* CD 



5/ 



v 7 a* 2 ax 



Km toj 



(22) 



20 Modifying the integration such that the integration is non-dimensionalized over the 
length L: 
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pA 0 



{x)\ I 0 pi 



4c 



^ 

V ; di ~dx 



rf+2mfl 0 
dx 



5/ 



(23) 



The length L cancels out, yielding non-dimensional modal mass, damping, Coriolis, 
5 and stiffness matrices mr, dr, cr, kr. 



mr„ 





~pA 0 "I ' 


K(*)j 


0 pi] 



*,fe)J 



(24) 
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(4-4) 



^0 5 £l 



[-Hi- 



(25) 



£//• „ + C7\„ „ = L / < 
m,n m,n i x y L j Vp 



rf + 2mfi 0 

o o 



(26) 
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(3 



(27) 



The modal differential equation of motion for flow can be reduced to matrix 



form: 
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[>ZZZ:]{^ (28) 



The modal response \ti(t_)\ may be considered to be harmonic, such that equation 
(28) is modified to equation (29) (for a beam with fixed ends, it would be preferable 
5 to drive the sensor conduit at co = co } ): 



(-6) 2 [mr] + 7^[^ + cr] + [k])fc(0} = {^} (29) 

10 The quantity in the parentheses on the left may be considered a "modal dynamic 
stiffness" analogous to that described above with reference to equations (6) and 
(7). A non-dimensional modal response can be calculated by taking the inverse of 
the modal dynamic stiffness (the modal frequency response function) and bringing 
it over to the right hand side to yield: 

15 

{^(/)} = (-^[ ZZ2 r] + jg)[dr^cr\^kr\f {N} (30) 

At this point, a simplifying assumption may be made that the beam acts like a 
Euler/Bemoulli beam, i.e. p£^A^S^0. For a Euler/Bernoulli beam, the modal 

20 stiffness may be expressed as: 



r*3 2 O m (x) d 2 O r (x) , 

] a - — ft r V Ljhr 

°" M 'K dx 2 dx 2 

. rV L d 2 o^(x)d 2 ^(x) 

t^ms ~" jxj y ' ^ aX 



(31) 



S Vl dx 2 



This allows the non-dimensional constants to be pulled out from the integration 
25 such that equation (31) simplifies to equation (32), where the a indicates a 
constant matrix due to integration of the mode shapes: 



fete)} = (~£ 2 pAm\+ ;^(d[^-]+2z[£r]) +[£:]) ' {K} (32) 
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A conduit having substantially fixed ends can be assumed to have a 
response that approximates that of a beam with fixed ends at a first bend mode 
frequency. Therefore, the physical response of the conduit can be approximated 
5 as if it were the sum of the product of fixed beam mode shapes times the fixed 
beam modal response, which means that mode shapes may be assumed to be 
effectively held constant at the excitation frequency. Therefore the modal mass, 
damping, Coriolis and stiffness matrices [mc],[A:],[s:]» can be assumed to 

be constant. Assuming the modal excitation doesn't change, the only variables in 
10 equation (32) that change are pA, d, and m . It may also be assumed that gy = 1 , 
as the first bending mode is typically driven in a straight-tube Coriolis meter. 
Theoretically, mass flow can be measured in modal domain because {n(t)} is 

complex with flow: 



PA = (( A4 + P^f ) L ) °>* JT 
15 , (33) 



25 



L 3 



m = 



_ Pf A f v f 



EI, 



where (p,A, + p f A f )L ^ bfse ° 2 ^ represents the mass of the flow tube and the fluid 
inside it, co x is the drive frequency (assuming driving of the first bending mode), and 
^L-lM}^ i S proportional to the stiffness of the beam. 

20 Expanding the non-dimensional parameters and factoring out the common 

variables yields: 



M(34) 
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If fluid density stays substantially constant in equation (34) with relation to the total 

F J / 

mass per unit length, the calibration bias is proportional to yl 3 . If this 



assumption is made, the remaining variables in equation (34) are EI t , P f A f v f , cc, 
and C , where co is the excitation frequency, p f A f v f is the mass flow rate, and 
EJ t is the flexural rigidity of the conduit. In a flow sensor application, the length 
typically is known and fixed by design. If the flow sensor conduit is constrained to 
substantially fixed boundary conditions, the mode shapes and modal frequencies 
may also be known. Therefore, the modal mass, damping, Coriolis and stiffness 
matrices [mr] . [</r] , [cr] , and [fa:] can be assumed to be constant. So the physical 
response can be transformed to any number of physical response points, p , where 
r is the number of modes. As long as the conduit ends are substantially fixed, the 
modal matrix used to transform the modal response back to a physical response is 
known: 



The physical response is complex. A mass flow measurement may be made from 
the physical response using time difference or phase measurements as described 
above. 

From equation (34) and (35), the only unknowns are the mass/unit length 
and El t . If a simplifying assumption that the sensor conduit is an Euler beam is 
made, the eigenvalues, X r , of the mode shapes may be determined. Therefore, if 
El t is measured, a calibration factor may be determined to a substantially high 
degree of accuracy, for example, <1% error. For an Euler/Bernoulli beam with 
fixed/fixed boundary conditions, the eigenvalues ( X r ) may be known. The 

pA 

frequencies of at least the drive mode are typically known. Therefore — can be 
estimated from the frequency and eigenvalue information as follows: 






(36) 
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For a more accurate sensor, secondary influences of shear stiffness and 
rotary inertia that effect the flow conduit's frequencies and boundary conditions may 
be considered. Tension effects may also be considered. The following shows how 
to back out the first four eigenvalue parameters in equation (20). Equation (20) 
provides the following non-dimensional relationships: 



pA 

— El 



kGAL 2 



EI, 



— EI, 



EI, 



(37) 



1 0 From equation (37), an equation for the frequencies of the modes can be derived, 
using the real normal mode assumptions: 



det 



-(S + B)^ B-X^-pho 2 



= 0 



(38) 



15 Equation (38) can be rearranged in terms of the eigenvalues, modal frequencies, 
and the non-dimensional eigenvalue parameters, yielding: 



20 



M r 2 -I2, 2 =[l o r 2 } 



2(A + S)^ + plj-4pA 

-(d+s) ' 



= [l a> 2 a>;]{Vlf 



= [l v r 2 ]{V2} 
(39) 



For a fixed/fixed beam, or a beam with other quantifiable boundary conditions, the 
left-hand side of equation (39) may be known and considered constant for many 
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modes. If modal frequencies for at least three modes are measured, one can solve 
for pA, S, pi , and B : 



S = -V2 0 = -JvT 0 

, V\,+2V2 n V2 
pA = 



0' ^1 



— 4 



V2 

B = -2pA 2 



2pA 



— V2 2 -V\ 3 — 



K V2 2 d -VlJ V2 2 '-V1 3 



pI = -V2 2 - 



pA 



(40) 



10 



The eigenvalue parameters may be broken down into pA,pI,EI„kGA, , and 
S using a modal mass determined from frequency response function (FRF) 
measurements. One way to measure pA and pi is to measure the modal mass of 
the flow conduit. A circle fit method from a limited amount of FRF measurements 
can provide an estimate of the mode's eigenvalue and residue. The modal mass 
m r of an r-th mode may be calculated from the residue R gqr , where qqr refers to a 

driving point response of the r-th mode: 



15 



R =a , 2 _Z^_ 
qqr qv 2m r o> r 

2 —j 



(41) 



m, = O 



qv 2Ra 

qqr r 



A least squares estimate of pA and pi may then be generated: 



V m ^l etmral = P A [^hpl[An r ) (42, 

20 Now [rhr 0 ] and [mr v ] can be considered as near constant (but not diagonal) 

matrices, and may be pulled from equation (42). A least squares estimate for pA 
and pi may be made by stacking the column vectors and taking the pseudo 
inverse: 
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(43) 



measured 



I, 



Generally, pA and pi will change as the density of the fluid changes. At the time 
5 of installation, the fluid density will probably be unknown. However the desired 
parameters EI t and kGA t are generally time invariant, although E and G change 
slightly with temperature and l t and A t may change with pressure. Therefore, at the 
time of installation, pA and pi can be measured and used to calculate EI t and 

kGA t . 

10 With this information in hand, an in situ estimate of a calibration factor may 

be made. pA may change the calibration bias if mass flow measurements are 
being made from phase measurements. Such a fluid density effect on calibration 
bias can be reduced by using time difference measurements instead of phase 
measurements. 

15 A potential disadvantage of the above technique is that measurement of the 

modal mass generally requires determination of a driving point response, e.g., a 
response at the location at which the flow conduit is driven. Such a driving point 
response is not required, however, to measure flow. Therefore, if it is undesirable 
to measure driving point response during operation of a flow sensor, one could use 

20 a motion transducer at the driving point for installation/calibration purposes and 
then remove the transducer (or ignore its output) during normal operation. 
Exemplary Calibration Apparatus and Operations 

FIG. 3 illustrates an exemplary implementation of a calibration factor 
generator circuit 300 according to embodiments of the present invention. The 

25 calibration factor generator circuit 300 is operative to determine a calibration factor 
325 for a flow sensor. As shown, the calibration factor generator circuit 300 
includes a normal modal dynamic characterizer circuit 310 that is operative to 
generate a normal modal dynamic characterization of conduit motion, e.g., a 
representation of a differential equation of motion in a modal domain comprising a 

30 plurality of normal modes. The calibration factor generator circuit 300 also includes 
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a calibration factor determiner circuit 320 that is operative to determine a calibration 
factor 325 from the normal modal dynamic characterization 315. For example, the 
calibration factor determiner circuit 320 may be operative to generate a solution to 
the modal differential equation of motion from the normal modal characterization 
5 315, and to determine the calibration factor from this solution, e.g., by generating a 
time difference or phase relationship estimate from the solution, and determining 
the calibration factor 325 therefrom. The calibration factor generator circuit 300 
further includes an interface circuit 330 operative to convey the determined 
calibration factor 325 to a mass flow estimator 360 of a mass flow sensor. 

10 As illustrated in FIG. 3, the normal modal dynamic characterizer circuit 310 

may generate the normal modal characterization 315 using a mode shape function 
O generated by a mode shape function generator circuit 350. The mode shape 
function generator circuit 350 generates the mode shape function O from motion 
signals 301 that represent motion of the conduit of the mass flow sensor, as 

15 described in the aforementioned United States Patent Application Serial Number 
09/941,332, the normal modal dynamic characterizer circuit 310 may also use 
modal mass, damping and stiffness terms [mr], [dr], [kr] generated by a modal 
analyzer circuit 340, or may generate such information using the mode shape 
function O , as described above with reference to equations (24), (25), and (27). In 

20 performing modal analysis of the flow sensor conduit, the modal analyzer circuit 
340 may, for example, generate drive signals 303 to excite the mass flow sensor 
conduit and may process motion signals 302 representing motion in response to 
this excitation using, for example, conventional modal analysis techniques. The 
normal modal dynamic characterizer circuit 310 may also use a spatial Coriolis 

25 matrix [C] to generate a modal Coriolis matrix term of the normal modal dynamic 

characterization 315, as described with reference to equation (3), or may determine 
such a modal Coriolis term directly from the mode shape function <D . 

As conceptually illustrated in FIG. 3, the calibration factor generator circuit 
300 may comprise a separate unit, or may be combined with the modal analyzer 

30 circuit 340 and/or the mode shape function generator circuit 350 and/or the mass 

flow estimator circuit 360. For example, all or part of the calibration factor 

generator circuit 300 may be integrated with the mass flow estimator circuit 360 are 

integrated in a common data processor, such as a microcontroller, microprocessor 

or digital signal processor (DSP) that forms part of a mass flow sensor apparatus, 

30 
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such that, for example, the normal modal dynamic characterizer circuit 310, the 
calibration factor determiner circuit 320, the interface circuit 330 and the mass flow 
estimator circuit 360 are embodied as software objects, modules or the like that are 
configured to cooperatively execute on the common data processor. Alternatively, 
5 all or part of the calibration factor generator circuit 300 may be integrated with the 
mode shape function generator circuit 350 and/or the modal analyzer circuit 340 in, 
for example, a common data processing apparatus used in a piece of test or 
process control equipment, such that, for example, these circuits are embodied as 
cooperatively executing objects, modules or the like. It will be appreciated that, in 

10 general, the calibration factor generator circuit 300 may be implemented using 

special purpose hardware, software or firmware executing on a general or special- 
purpose data processor, or combinations thereof. 

FIGs. 4-6 are flowchart illustrations of exemplary operations according to 
various embodiments of the present invention. Those skilled in the art will 

1 5 understand that the operations of these flowchart illustrations may be implemented 
using computer instructions. These instructions may be executed on a computer or 
other data processing apparatus, such as may be used to implement the calibration 
factor generator circuit 300 of FIG. 3, to create an apparatus (system) operative to 
perform the illustrated operations. The computer instructions may also be stored as 

20 computer readable program code on a computer readable medium, for example, an 
integrated circuit memory, a magnetic disk, a tape or the like, that can direct a 
computer or other data processing apparatus to perform the illustrated operations, 
thus providing means for performing the illustrated operations. The computer 
readable program code may also be executed on a computer or other data- 

25 processing apparatus to cause the apparatus to perform a computer-implemented 
process. Accordingly, FIGs. 4-6 support apparatus (systems), computer program 
products and methods for performing the operations illustrated therein. 

Those skilled in the art will appreciate that the present invention may be 
implemented a number of other ways than the embodiments described herein. For 

30 example, computations described herein may be implemented as separate 

computations, or may be combined into one or more computations that achieve 
equivalent results. The functions described herein may, in general, be 
implemented using digital and/or analog signal processing techniques. Those 
skilled in the art will also appreciate that, although the present invention may be 
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embodied within an apparatus such as a Corioiis mass flowmeter, or as methods 
which may be performed by such apparatus, the present invention may also be 
embodied in a apparatus configured to operate in association with a flowmeter or 
sensor apparatus, such as in process control apparatus. It will also be appreciated 
5 that the present invention may be embodied in an article of manufacture in the form 
of computer-readable instructions or program code embodied in a computer 
readable storage medium such as a magnetic disk, integrated circuit memory 
device, magnetic tape, bubble memory or the like. Such computer program code 
may be executed by a computer or other data processor and executed responsive 

10 to motion signals supplied from motion transducers operatively associated with a 
flow sensor conduit or similar structure. 

FIG. 4 illustrates exemplary operations 400 for generating a calibration factor 
for a flow sensor according to embodiments of the present invention. A mode, 
shape function, e.g., such as the mode shape functions described above with 

15 reference to equation (2), are determined from measured response data, from 
motion signals such as those produced by the motion transducers of FIG., 1 (Block 
410). A normal modal dynamic characterization is then generated from the mode 
shape function (Block 420). For example, a modal differential equation of motion 
including a modal mass term, a modal stiffness term, a modal damping term and a 

20 modal Corioiis term may be generated by determining modal mass and stiffness 
matrices using modal analysis techniques and generating a modal Corioiis matrix 
from a discrete physical Corioiis matrix for a predetermined mass flow, as 
described above with reference to equations (3) - (15), which may be viewed as a 
determination of orthogonality of discrete mode shapes with respect to the spatial 

25 Corioiis characterization (e.g., the physical Corioiis matrix). Alternatively, a modal 
Corioiis matrix may be determined by determining an orthogonality of a continuous 
mode shape function (e.g., with respect to itself and other mode shape functions) 
for a predetermined mass flow rate, as described above with reference to equation 
(26), and modal mass and stiffness matrices may be determined by a similar 

30 orthogonalization, e.g., equations (24), (25), and (27), or by conventional modal 
analysis techniques. A calibration factor may then be determined from the normal 
modal characterization (Block 430). For example, a modal differential equation of 
motion including modal mass, stiffness, damping and Corioiis terms may be solved 
and translated to a spatial domain to yield an estimated spatial response from 
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which a calibration factor may be determined, as described above with reference to 
equations (3), (6) - (8) and (32) - (35). 

FIG. 5. illustrates exemplary operations 500 for determining a flow sensor 
calibration factor according to some embodiments of the present invention. Modal 
5 mass and stiffness terms of a modal differential equation of motion, e.g., the [mr] 
and [/cr] matrices of equation (6), are determined using, for example, using 
conventional modal analysis techniques (Block 510). A modal Coriolis term of the 
equation of motion, e.g., the [cr] matrix of equation (6), is then determined from a 
spatial Coriolis characterization and a mode shape function, e.g., from a physical 

10 Coriolis matrix [C] as described in equation (3) (Block 520). A calibration factor is 
then determined from a normal modal dynamic characterization of conduit motion 
as a function of the modal mass, stiffness and Coriolis terms, e.g., from the modal 
differential equation of motion (Block 530). 

FIG. 6 illustrates exemplary operations 600 for determining a flow sensor 

1 5 calibration factor according to other embodiments of the present invention. Modal 
mass and stiffness terms of a modal differential equation of motion, e.g., the [mr] 
and [k}-\ matrices of equation (34), are determined using, for example, conventional 
modal analysis techniques or by determining orthogonality of a mode shape 
function as described with reference to equations (24) and (25) (Block 610). A 

20 modal Coriolis term of the equation of motion, e.g., the [cr] matrix of equation (34), 
is then determined by determining orthogonality of a mode shape function for a 
predetermined mass flow rate, e.g., as described with reference to equation (26) 
(Block 620). A calibration factor is then determined from a normal modal dynamic 
characterization of conduit motion as a function of the modal mass, stiffness and 

25 Coriolis terms, e.g., from a solution of the modal differential equation of motion 
(Block 630). 

As noted above, simplifying assumptions as to the modal behavior of a flow 
sensor conduit may be made if the motion of the conduit is properly constrained to 
provide predetermined boundary conditions. In particular, substantially fixed, 
30 homogeneous boundary conditions for at least one mode of the conduit may be 
achieved by affixing the ends of a sensor conduit. 

For example, as shown in FIG. 7, a flow sensor 700 may comprise a conduit 
703 that is driven (vibrated) by an actuator 706. Motion sensors 705 generate 
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motion signals 707 that represent motion of the conduit 703 as a material 708 flows 
therethrough . A mass flow determiner circuit 710 generates a mass flow estimate 
715 from the motion signals 707 using a calibration factor 755 derived from a 
normal modal dynamic characterization of the sensor 700. 
5 In particular, as described above, the calibration factor 755 may be derived 

assuming a predetermined boundary condition for a normal mode of the conduit 
703, for example, assuming a fixed/fixed boundary condition for a first bending 
mode driven by the actuator 706. As conceptually illustrated in FIG. 7, this 
boundary condition may be approximated by fixing locations 703A, 703B of the 

10 conduit 703 to a substantially rigid structure 730. It will be appreciated that this 
fixation of the locations 703A, 703B may be achieved in number of different ways. 
For example, the locations 703A, 703B of the conduit 704 may be affixed to a rigid 
structure using clamps, welds or other fastening techniques. 

It will be appreciated that according to embodiments of the present invention, 

15 such an approach can facilitate implementation of a simplified straight tube^Goriolis 
flowmeter that does not require a mechanically-tuned balanced beam as used in 
many conventional straight tube sensors. It will be further appreciated that the 
above-described techniques can also enable practical implementation of a "clamp- 
on" flow sensor, I.e., a flow sensor implemented by transducers, actuators and 

20 associated circuitry with a length of pipe or other conduit that is part of a processing 
apparatus, e.g., an existing pipeline in a refinery, processing plant, or other 
installation. 

An example of a straight tube flow sensor configuration according to 
embodiments of the present invention is illustrated in FIG. 8. In FIG. 8, a flow 

25 sensor 800 includes a conduit 803 surrounded by and attached to a case 804. The 
conduit 803 is configured to receive and discharge material from send to, for 
example, a pipeline (not shown) attached at flanges 802. The sensor 800 further 
includes one or more actuators 806 that induce relative motion between the conduit 
803 and the case 804. A plurality of motion transducers 805 generate motion 

30 signals 807 that represent relative motion of the conduit 803 with respect to the 
case 804. 

The sensor interface circuit 810 includes a mass flow determiner circuit 812 
that is operative to receive the motion signals 807 and to generate a mass flow 
estimate 813 therefrom according to a calibration factor 855. The sensor interface 
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circuit 810 further includes a driver circuit 814 operative to control the one or more 
actuators 806. 

The calibration factor 855 is generated by a normal modal dynamic 
characterizing calibration circuit 850. The calibration circuit 850 derives the 
5 calibration factor 855 from a normal modal dynamic characterization of motion of 
the conduit 803. As shown by the dashed line connections, the calibration circuit 
850 may be operative to receive motion signals 807 generated by the transducers 
805 and/or to generate drive signals to drive the one or more transducers 806. 
This capability may be used in determining mode shape functions and terms of a 

1 0 modal differential equation of motion, and other parameters needed for 
determination of the calibration factor 855. 

As described above, the calibration factor 855 may be derived based on an 
assumption that the conduit 803 will be constrained to a predetermined boundary 
condition, such as a fixed end boundary condition for a first bending mode of the 

1 5 conduit 803. To approximate such a predetermined boundary condition, the motion 
of the conduit 803 is constrained. For example, as shown conceptually in FIG. 8, 
motion near the flanges 802 may be constrained to a near fixed condition by 
affixing the conduit 803 to a rigid structure such that motion near the flanges 802 is 
substantially reduced. It will be appreciated that such fixation may be achieved in a 

20 number of different ways. For example, motion near the flanges 802 may be 

constrained by clamping, welding or otherwise fixing a pipeline connected to the 
flanges 802 to a rigid structure. Fixation may also be achieved, for example, by 
affixing the case 804 to a rigid structure, assuming the case 804 is sufficiently rigid 
to constrain motion near the flanges 802. 

25 It will be further appreciated that the sensor interface circuit 810 may 

implemented in a number of different ways. For example, the sensor interface 
circuit 810 may be implemented in electronics packages with the case 804. 
Alternatively, all or portions of the sensor interface circuit 810 may be positioned 
away from the case 804 and components therein. For example, example, all or 

30 portions of the sensor interface circuit 810 may be included in remote process 
control equipment connected to the one or more actuators 806 and/or the motion 
transducers 805. 

It will be further appreciated that the calibration circuit 850 may be 
implemented in a number of different ways. For example, all or portions of the 
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calibration circuit 850 may be integrated with the sensor interface circuit 810. All or 
portions of the calibration circuit 850 may also be included, for example, in remotely 
located process control or other equipment designed to interface with the sensor 
800 in a plant or other environment. All or portions of the calibration circuit 850 
5 may also be provided, for example, in test equipment (field or factory) that is 

configured to be connected to the sensor 800 for calibration purposes, but removed 
or deactivated for normal operation of the sensor 800. The sensor interface circuit 
810 may also share components, such as processing circuits and driver circuits, 
with the calibration circuit 850. 

10 FIG. 9 illustrates an example of a "clamp-on" flow sensor configuration 

according to other embodiments of the present invention. A pipeline 903 is 
clamped to a rigid structure 930 at first and second locations 903A, 903B by first 
and second pipeline clamps 911 A, 91 1B. A plurality of motion transducers 905, 
e.g., inertial or relative transducers, are configured to engage the pipeline 903 at 

15 respective locations, and produce motion signals 907 representative of motion of 
the conduit 903. A sensor interface circuit 910 includes a mass flow determiner 
circuit 912 that is operative to process the motion signals 907 according to 
calibration factor 955 derived from a normal modal dynamic characterization of 
motion of the pipeline 903. One or more actuators 906 are also configured to 

20 engage the pipeline 903, and are driven by a driver circuit 914 included in the 
sensor interface circuit 910. 

The calibration factor 955 is generated by a normal modal dynamic 
characterizing calibration circuit 950. The calibration circuit 950 derives the 
calibration factor 955 from a normal modal dynamic characterization of motion of 

25 the conduit 903. As shown by the dashed line connections, the calibration circuit 
950 may be operative to receive motion signals 907 generated by the transducers 
905 and/or to generate drive signals to drive the one or more transducers 906. 
This capability may be used in determining, for example, mode shape functions and 
terms of a modal differential equation of motion, and other parameters needed for 

30 determination of the calibration factor 955. 

The calibration factor 955 may be derived based on an assumption that the 

conduit 903 will be constrained to a predetermined boundary condition, such as a 

fixed end boundary condition for a first bending mode of the conduit 903. To 

approximate such a predetermined boundary condition, the motion of the conduit 

36 
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903 is constrained. For example, as shown in FIG. 8, motion at the first and 
second spaced apart locations 903A, 903B of the conduit 903 may be constrained 
by the clamps 91 1 A, 91 1 B that attach the conduit 903 to the rigid structure 930. It 
will be appreciated that such fixation may be achieved in other ways. 
5 The sensor interface circuit 910 may implemented in a number of different 

ways. For example, the sensor interface circuit 910 may be implemented in an 
electronics package configured to be positioned near the actuators 905 and the 
driver 906. Alternatively, all or portions of the sensor interface circuit 910 may be 
positioned away from the transducers 905 and the one or more actuators 906. For 

10 example, example, all or portions of the sensor interface circuit 910 may be 
included in remote process control equipment connected to the one or more 
actuators 906 and/or the motion transducers 905 via communications links. 

It will be further appreciated that the calibration circuit 950 may be 
implemented in a number of different ways. For example, all or portions of the 

15 calibration circuit 950 may be integrated with the sensor interface circuit 910. All or 
portions of the calibration circuit 950 may also be included, for example, in process 
control or other equipment All or portions of the calibration circuit 950 may also be 
provided, for example, in test equipment (field or factory) that is connected to the 
actuators 906 and transducers 905 for calibration purposes, but removed or 

20 deactivated for normal operation. The sensor interface circuit 910 may also share 
components, such as signal processing circuits and driver circuits, with the 
calibration circuit 950. 

FIG. 10 illustrates another example of a clamp-on flow sensor configuration 
according to embodiments of the present invention. A pipeline 1003 is clamped to 

25 a rigid structure 1030 at first and second locations 1003A, 1003B by first and 
second pipeline clamps 10I1A, 1011B. The clamps 1011A, 1011B are also 
attached to a rigid structure 1004, such as a rigid beam. A plurality of relative 
motion transducers 1005 are configured to engage the pipeline 1003 at respective 
locations, and provide motion signals 1007 representative of motion of the conduit 

30 1003 relative to the structure 1004. A sensor interface circuit 1010 includes a mass 

flow determiner circuit 1012 that is operative to process the motion signals 1007 

according to calibration factor 1055 derived from a normal modal dynamic 

characterization of motion of the pipeline 1003. One or more actuators 1006 are 

also configured to engage the pipeline 1003, and to drive the conduit 1003 relative 
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to the structure 1004 responsive to a driver circuit 1014 included in the sensor 
interface circuit 1010. 

The calibration factor 1055 is generated by a normal modal dynamic 
characterizing calibration circuit 1050. The calibration circuit 1050 derives the 
5 calibration factor 1055 from a normal modal dynamic characterization of motion of 
the conduit 1003. As shown by the dashed line connections, the calibration circuit 
1050 may be operative to receive motion signals 1007 generated by the 
transducers 1005 and/or to generate drive signals to drive the one or more 
transducers 1006. This capability may be used, for example, in determining mode 

10 shape functions and terms of a modal differential equation of motion, and other 
parameters needed for determination of the calibration factor 1055. 

The calibration factor 1055 may be derived based on an assumption that the 
conduit 1003 will be constrained to a predetermined boundary condition, such as a 
fixed end boundary condition for a first bending mode of the conduit 1003. To 

15 approximate such a predetermined boundary condition, the motion of the conduit 
1003 is constrained. For example, as shown conceptually in FIG. 8, motion at first 
and second spaced apart locations 1003A, 1003B of the conduit 1003 may be 
constrained by the clamps 1011 A, 101 1B and beam 1004 are attached to a rigid 
structure 1030. It will be appreciated that such fixation may be achieved in other 

20 ways. 

The sensor interface circuit 1010 may be implemented in a number of 
different ways. For example, the sensor interface circuit 1010 may be implemented 
in an electronics package configured to be positioned near the actuators 1005 and 
the driver 1006. Alternatively, all or portions of the sensor interface circuit 1010 

25 may be positioned away from the transducers 1005 and the one or more actuators 
1006. For example, all or portions of the sensor interface circuit 1010 may be 
included in remote process control equipment connected to the one or more 
actuators 1006 and/or the motion transducers 1005 via communications links. 
It will be further appreciated that the calibration circuit 1050 may be 

30 implemented in a number of different ways. For example, all or portions of the 

calibration circuit 1050 may be integrated with the sensor interface circuit 1010. All 

or portions of the calibration circuit 1050 may be included, for example, in process 

control or other equipment. All or portions of the calibration circuit 1050 may also 

be provided, for example, in test equipment (field or factory) that is connected to the 
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actuators 1006 and transducers 1005 for calibration purposes, but removed or 
deactivated for normal operation. The sensor interface circuit 1010 may also share 
components, such as signal processing circuits and driver circuits, with the 
calibration circuit 1050. 
5 According to other embodiments of the present invention, predetermined 

boundary conditions for a flow sensor conduit may also be approximated by using 
mode shape control techniques that constrain motion of a conduit of a mass flow 
sensor. Such mode shape control techniques are described, for example, in United 
States Patent Application Serial Number 09/942,189, entitled "Sensor Apparatus, 

10 Methods and Computer Program Products Employing Vibrational Shape Control," 
to Wheeler, filed concurrently herewith and incorporated herein by reference in its 
entirety. In particular, the aforementioned application describes how actuators in 
operative association with a flow sensor conduit may be used to approximate a 
fixed boundary condition for at least one vibrational mode of the conduit, for 

1 5 example, a first bending mode for a straight-tube type conduit sensor, 

FIG. 11 illustrates a mass flow sensor apparatus 1 100 using such mode 
shape control according to embodiments of the present invention. The apparatus 
1 100 includes a conduit 1 103. A plurality of actuators 1 106 are operative to apply 
force to the conduit at a plurality of locations. A plurality of motion transducers 

20 1 105 are also operatively associated with the conduit 1 103 and generate motion 
signals 1107 responsive to motion of the conduit 1103. 

Two or more of the motion signals 1107 may be processed by a mass flow 
determiner circuit 1 1 12 of a sensor interface circuit 1 1 10 to generate a mass flow 
estimate 1113 using a calibration factor 1 155 derived from a normal modal dynamic 

25 characterization of motion of the conduit 1 103 based on an assumption of a 

predetermined boundary condition for a mode (e.g., a first bending mode) of the 
conduit 1103. For example, the assumed predetermined boundary condition may 
comprise a substantial lack of motion (e.g., in an inertial frame) at the locations 
1103A, 1103Bofthe conduit 1103. 

30 The sensor interface circuit 1110 further includes a mode shape control 

circuit 1 1 14 that is operative to receive some or all of the motion signals 1 107 and 
to generate the drive signals 1 109 responsive thereto. For example, the mode 
shape control circuit 1120 may be operative to vibrate the conduit 1103 in a first 
bending while constraining motion of the conduit 1103 such that the predetermined 
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boundary condition upon which the calibration factor 1155 is based is 
approximated. For example.the mode shape control circuit 1 120 may drive the 
actuators 1 106 such that a first bending mode of the conduit 1 103 is driven while 
maintaining substantially zero motion in the first bending mode at the conduit 
5 locations 1 1 03A, 1 1 03B. 

The calibration factor 1155 is generated by a normal modal dynamic 
characterizing calibration circuit 1 150. The calibration circuit 1 150 may be 
operative to receive motion signals 1 107 generated by the transducers 1 105 and/or 
to generate drive signals to drive the one or more transducers 1 106 (connections 

1 0 are not shown in FIG. 1 1 for purposes of clarity of explanation). This capability may 
be used, for example, in determining mode shape functions and terms of a modal 
differential equation of motion, and other parameters needed for determination of 
the calibration factor 1 155. 

The sensor interface circuit 1110 may be implemented in a number of 

1 5 different ways. For example, the sensor interface circuit 1110 may be implemented 
in an electronics package configured to be positioned near the actuators 1105 and 
the drivers 1 106. Alternatively, all or portions of the sensor interface circuit 1110 
may be positioned away from the transducers 1 105 and the one or more actuators 
1 1 06. For example, example, all or portions of the sensor interface circuit 1110 

20 may be included in remote process control equipment connected to the actuators 
1006 and/or the motion transducers 1 1 05 via communications links. 

It will be further appreciated that the calibration circuit 1 150 may be 
implemented in a number of different ways. For example, all or portions of the 
calibration circuit 1 150 may be integrated with the sensor interface circuit 1110. All 

25 or portions of the calibration circuit 1 1 50 may be included, for example, in process 
control or other equipment located distant from the conduit 1 1 03. All or portions of 
the calibration circuit 1 150 may also be provided, for example, in test equipment 
(field or factory) that is connected to the actuators 1 106 and transducers 1 105 for 
calibration purposes, but removed or deactivated for normal operation. The sensor 

30 interface circuit 1110 may also share components, such as signal processing 
circuits and driver circuits, with the calibration circuit 1 150. 

It will be further understood that, similar to the configurations described 
above with reference to FIGs. 8-10, the active vibration control approach described 
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above with reference to FIG. 1 1 may be used to implement straight-tube flow 
sensors, including integrated and clamp-on configurations. 

FIG. 12 illustrates exemplary operations 1200 according to yet other aspects 
of the present invention. A calibration factor for a mass flow sensor is determined 
from a normal modal dynamic characterization of the sensor's conduit assuming a 
predetermined boundary condition, such as a fixed end boundary condition for one 
or more normal modes of the conduit (Block 1210). Motion signals representing 
motion of the conduit are generated while constraining motion of the conduit to 
approximate the predetermined boundary condition, e.g., by attaching appropriate 
locations of the conduit to a rigid structure and/or by using shape control (Block 
1220). The motion signals so generated are processed according to the calibration 
factor to generate a mass flow estimate (Block 1 230). 

In the drawings and specification, there have been disclosed typical 
embodiments of the invention and, although specific terms are employed, they are 
used in a generic and descriptive sense only and not for purposes of limitation, the 
scope of the invention being set forth in the following claims. 
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THAT WHICH IS CLAIMED IS: 

1 . A method of determining mass flow in a conduit configured to contain 
a material, the method characterized by the steps of: 

determining a calibration factor from a normal modal dynamic 
characterization of motion of the conduit; 

generating a plurality of motion signals representing motion of the conduit at 
a plurality of locations on the conduit; and 

processing the motion signals according to the determined calibration factor 
to produce a mass flow estimate. 



2. A method according to Claim 1 , wherein determining a calibration 
factor comprises: 

generating a solution to a modal differential equation of motion in a modal 
domain comprising a plurality of normal modes of motion of the conduit; and 
1 5 generating the calibration factor from the generated solution of the modal 

differential equation of motion. 

3. A method according to Claim 1 , wherein the calibration factor relates 
mass flow rate to a spatio-temporal relationship among the plurality of motion 

20 signals. 



4. A method according to Claim 3, wherein the calibration factor relates 
mass flow rate to one of a time difference or a phase relationship. 

25 5. A method according to Claim 1 , wherein the normal modal dynamic 

characterization characterizes motion of the conduit as a function of a modal 
Coriolis term that describes coupling among a plurality of normal modes of the 
conduit responsive to mass flow in the conduit 

30 6. A method according to Claim 5, wherein determining a calibration 

factor is preceded by determining the modal Coriolis term from a mode shape 
function that describes motion of the conduit in a normal mode as a function of 
location on the conduit. 
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7. A method according to Claim 5, wherein determining a calibration 
factor is preceded by determining the modal Coriolis term from a spatial Coriolis 
characterization that describes motion of a plurality of discrete locations of the 
conduit in a spatial domain in response to a predetermined mass flow. 

5 

8. A method according to Claim 6, wherein the mode shape function is a 
function of an eigenvalue and a boundary condition. 

9. A method according to Claim 6, further comprising determining the 
1 0 mode shape function. 

10. A method according to Claim 9, wherein determining the mode shape 
function comprises determining the mode shape function from a predetermined 
eigenvalue and a predetermined boundary condition. 

15 

11. A method according to Claim 6, wherein determining a calibration 
factor is preceded by determining an orthogonality of a mode shape function for a 
predetermined mass flow in the conduit, the mode shape function describing 
motion of the conduit in a normal mode as a function of location on the conduit 

20 

12. A method according to Claim 1 1 , further comprising determining the 
mode shape function. 

13. A method according to Claim 1 , wherein determining a calibration 
25 factor comprises: 

generating an estimated spatial response for a plurality of locations of the 
conduit from the normal modal dynamic characterization; and 

generating the calibration factor from the estimated spatial response. 

30 14. A method according to Claim 13, wherein generating the calibration 

factor from the estimated spatial response comprises: 

determining a spatio-temporal relationship among movements at the plurality 
of locations; and 
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determining the calibration factorfrom the determined spatio-temporal 
relationship. 

15. A method according to Claim 1, wherein determining a calibration 

5 factor comprises generating the calibration factor from motion signals produced by 
the plurality of motion transducer 

16. A mass flow sensor comprising a conduit (103) configured to contain 
a material (108), a plurality of motion transducers (105) operative to generate 

10 motion signals representing motion of the conduit, at least one actuator (106) 

operative to excite the conduit and a mass flow estimator circuit (360) operative to 
estimate mass flow from the motion signals according to a calibration factor, the 
apparatus characterized by: 

a calibration factor generator circuit (300) operative to generate a calibration 

1 5 factor from a normal modal dynamic characterization of motion of the conduit, the 
calibration factor generator circuit including an interface circuit (330) operative to 
convey the generated calibration factor (325) to the mass flow estimator circuit. 

17. A mass flow sensor according to Claim 16, wherein the calibration 
20 factor generator circuit comprises a normal modal dynamic characterizer circuit 

(310) operative to generate a normal modal dynamic characterization of motion of 
the conduit; and 

18. A mass flow sensor according to Claim 16, wherein the normal modal 
25 dynamic characterization characterizes motion of the conduit as a function of a 

- modal Coriolis term that describes coupling among a plurality of normal modes of 
the conduit responsive to mass flow in the conduit. 

19. A mass flow sensor according to Claim 18, wherein the calibration 
30 factor generator circuit is operative to determine the modal Coriolis term from a 

mode shape function that describes motion of the conduit in a normal mode as a 
function of location on the conduit. 
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20. A mass flow sensor according to Claim 19, wherein the calibration 
factor generator circuit is operative to determine the mode shape function from a 
predetermined eigenvalue and a predetermined boundary condition. 

5 21 . A mass flow sensor according to Claim 16, wherein the calibration 

factor generator circuit is operative to generate an estimated spatial response for a 
plurality of locations of the conduit from the normal modal dynamic characterization 
and to generate the calibration factor from the estimated spatial response. 

10 22. A mass flow sensor according to Claim 16, wherein the calibration 

factor generator circuit is configured to receive motion signals from the plurality of 
motion transducers of the parameter sensor and operative to generate the 
calibration factor from the received motion signals. 

15 23. A mass flow sensor according to claim 16, wherein the mass flow 

estimator circuit is configured to receive the motion signals and operative to 
generate a mass flow estimate therefrom using a calibration factor derived from a 
normal modal characterization of motion of the conduit that assumes a 
predetermined boundary condition for the conduit and further comprising means 

20 (730A & 730B) for constraining motion of the conduit to approximate the 
predetermined boundary condition. 

24. A mass flow sensor according to Claim 23, wherein the mass flow 
estimator circuit is operative to determine one of a time difference or a phase 

25 relationship and to apply the calibration factor to the determined one of a time 
difference or a phase relationship to generate an estimate of mass flow rate: 

25. A mass flow sensor according to Claim 23, wherein the means for 
constraining motion of the conduit comprises means for applying a force to the 

30 conduit responsive to motion of the conduit to thereby induce a boundary condition 
that approximates the assumed boundary condition. 
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26. A mass flow sensor according to Claim 23 further comprising a shape 
control circuit configured to receive motion signals from at least two motion 
transducers and operative to control the at least one actuators responsively thereto. 
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